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We show that a continuous quantum non-demolition measurement of the energy of a nanome- 
chanical resonator can be achieved by monitoring the resonator with a quantum point contact via 
a Cooper-pair box. This technique can further be used to prepare highly non-classical states of two 
resonators, such as canonical phase-reference states, and so-called noon states. 
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Nanoscopic mechanical resonators can now be built 
with frequencies in the hundreds of Megahertz and qual- 
ity factors above 10 4 0, i, 0, i, i, i]. These resonators 
have potential applications in metrology Q and infor- 
mation processing [1] , and hold the promise of realizing 
quantum behavior in bulk mechanical devices for the first 
time. To fully exploit the quantum properties of nanores- 
onators one must be able to prepare them in nonclassical 
states, but generating the necessary nonlinear evolution 
is challenging. Our purpose here is twofold. First wc 
present a method for making a continuous QND (Quan- 
tum Non-Demolition) measurement of the energy of a 
nano-resonator. Such measurements are important be- 
cause they will prepare non-classical Fock states of the 
resonator, and allow the observation of quantum jumps. 
While these measurements have been a problem of con- 
siderable interest for a number of years, only two meth- 
ods to make such measurements in the solid-state exist 
to date [1, [H], and only the latter, suggested recently, 
appears feasible with current technology. (Wc note that 
two methods to monitor energy in a near-QND fashion 
have also been proposed recently ll|, [12J, as well as a 
promising QND scheme using an optical cavity [liij.) Sec- 
ond, we show that joint QND measurements of two res- 
onators can be used to prepare highly non-classical entan- 
gled states; in particular canonical phase-reference (CP) 
states [3 IH, till and so-called "noon" states [13, [11]. 
Both can be thought of as states of a single virtual oscil- 
lator whose number basis consists of the states of definite 
phonon number-difference between two oscillators. These 
states are very difficult to create via Hamiltonian evolu- 
tion, due to the high-order nonlinearities required [l9| . 
Both CP states and noon states have applications in high- 
precision phase measurement. The former also contain a 
large store of entanglement, and are therefore of potential 
use in information processing. Further, when combined 
with a universal procedure to create states of a single res- 
onator, such as that of Law and Eberly [201 ] . an energy 
measurement can be used to create any state of the vir- 
tual oscillator. While we will specifically consider nano- 



resonators in what follows, the methods we describe here 
can just as readily be applied to measurements and state- 
preparation of a superconducting striplinc resonator 21 1 . 
Methods to entangle solid-state qubits using measure- 
ments have also been proposed (see e.g. Refs. [22I l23j|). 

To perform a QND measurement on a nano-resonator 
one requires an interaction Hamiltonian with a probe sys- 
tem proportional to a^a where a is the annihilation op- 
erator for the resonator. Such an interaction can be ob- 
tained quite easily by placing a Cooper-pair Box (CPB) 
adjacent to the resonator and ensuring that the detun- 
ing between them, A, is large compared to their mutual 
interaction strength A [24j |. The resulting interaction is 
Hi u t = h/ia x a'a, where the charge basis for the CPB is 
the eigenbasis of a z , and p = A 2 /A. This kind of inter- 
action has already been used to obtain spectral measure- 
ments revealing the discrete energy levels of a supercon- 
ducting stripline resonator [2{| . 

To realize a continuous QND measurement we add one 
more element whose function is to perform a continuous 
measurement of the charge states of the CPB. This ele- 
ment can be either a quantum point-contact (QPC) [111 
or single-electron transistor (SET) placed adjacent to the 
CPB. (A QPC is more efficient than an SET, but as we 
discuss below efficiency is not essential for an effective 
QND measurement). Monitoring the current through ci- 
ther of these devices provides a continuous measurement 
of o z . The circuit diagram for the measurement scheme 
is shown in Fig. [1] The full stochastic master equation 
(SME) for the coupled resonator-CPB system, under the 



continuous measurement is 



dp = — i[H/h, p]dt — k[a Zl [a z , p]]dt 

+V2k[a z p + pa z - 2(a z )p](dr - (a z )dt), (1) 

where H = H[uiRa^ a + ioqo ~ z +uija x ] +Hi nt is the Hamilto- 
nian, is the resonator frequency, Eq is the CPB charg- 
ing energy, and wj is it Josephson tunneling frequency. 
The parameter k is the strength [28| of the measurement 
of a z and is determined by the coupling strength between 
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the SET/QPC and the CPB. This master equation is 
stochastic because it is driven by the continuous random 
stream of measurement results r(t). The increment of r 



in time dt is given by dr = (a z )dt + dW [27|, where dW 
is a Gaussian noise increment satisfying the Ito calculus 
relation 
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We show that this configuration realizes a continuous 
QND measurement of the energy of the resonator by sim- 
ulating the master equation Eq.(TT|). Starting the res- 
onator in a uniform superposition of the first ten number 
states, we find that the variance of the resonator's en- 
ergy decays to zero and remains there as required, with 
the resonator projected onto a single number state. In 
Fig. \2l&) we show the decay of the energy variance av- 
eraged over a thousand realizations. It is worth noting 
that the projection of the system onto an energy eigen- 
state is unaffected by any environmental disturbance to 
the CPB by virtue of the QND interaction. Decoherence 
can only affect the measurement by reducing the infor- 
mation extraction rate (by interfering with the rotation 
induced by the resonator). However, one expects this to 
be significant only if the decoherence rate is much greater 
than /i. 

The operating principle behind the continuous mea- 
surement is that the effective rotation rate u> x of the CPB 
about its x-axis depends directly on the number state of 
the resonator {u> x = u>j + KfMva). Because the SET- 
measurement is sensitive to the charge of the CPB (ct 2 ), 
it is also sensitive to the rate at which the z-cigenstate is 
changing, and hence the phonon number of the resonator. 
The measurement is QND by virtue of the CPB-resonator 
interaction. The projected phonon number can be deter- 
mined either from the full SET measurement record, or 
simply from the peak in its noise-spectrum (30l | . 

The rate at which information is extracted from the 
resonator, being the rate at which the resonator is pro- 
jected onto a phonon number state, is determined by the 
interaction strength with the CPB, /x, and the measure- 
ment strength k. For a given value of \i we now wish to 
determine the value of k that maximizes the projection 
rate. If we make k too small then we expect this to limit 
the rate of information extraction. On the other hand, 




FIG. 1: The circuit diagram for the QND measurement 
scheme: V and V g are applied voltages, and Ci and Ci are the 
capacitances connecting the three mesoscopic elements. The 
Cooper-pair box (indicated with a dashed box) is coupled to 
the nano-resonator and the single-electron transistor. 




FIG. 2: (a) The evolution of the phonon-number variance of 
the resonator for k — /i averaged over a thousand measure- 
ments. Inset: The average value of the phonon number for a 
single measurement, (c) This average variance as a function 
of k after measuring for a time t = 2-7r//x. 



making k too large will do the same: the information re- 
garding the energy is contained in the rate at which the 
CPB flips between z-states, and the quantum Zeno effect 
induced by the measurement will prevent this flipping. 
There will therefore be an optimal value of k for a given 
fi. We find this value by calculating the energy variance 
after a fixed measurement time, and averaging this over 
many simulations of the measurement. The results are 
presented in Fig. [DJb) , showing that the optimal choice 
is k as fj,. Realistic parameters are /i = 10 6 s _1 [TJ], and 
k can be much larger [3l[ . From Fig. [Ha) we see that 
this gives a measurement time of about 20fis. 



We now show that energy measurements are a powerful 
tool for preparing highly non-classical entangled states. 
In particular we show that they can be used to convert 
single-oscillator Schrodinger cat states, which are rela- 
tively easy to prepare, into two kinds of phase-reference 
states. This can be achieved with any detector capable 
of a joint-oscillator energy measurement. For our mea- 
surement scheme, we place two nano-resonators symmet- 
rically on cither side of the CPB. In this case the CPB 
interacts with both resonators simultaneously, and the in- 
teraction Hamiltonian is Hint = n<J x (a^ a + b^b), where b 
is the annihilation operator for the second oscillator. The 
measurement now extracts information about the sum of 
the energies of the two oscillators. The measurement will 
therefore project the two oscillators into a subspace with 
a fixed total phonon number. If the states of the two res- 
onators prior to the measurement are \ipi) = J2n c n\ n ) 
and |^2) = *52 n d n \n), respectively, then if the measure- 
ment projects onto the space with a total of N photons, 
the combined state of the resonators becomes 
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Each state in the sum, \n) at, is a state of definite phonon- 
number difference between the resonators. If each res- 



onator is prepared in the same (arbitrary!) coherent 
state prior to the measurement, then \ip) is an en- 
tangled state with the "N-bmomial" coefheients q n oc 
y/l/JjN — n)\n\). Fig[3Jb) shows a typical distribution 
\q n \ 2 for this case. 

We consider first preparing "noon" states of two oscil- 
lators, defined as \M) = \N)\0) + \0)\N) . These states are 
useful for very high precision measurements of a phase 
shift of one of the oscillators [32| (the precision scales 
linearly with N) . Consequently, there has been consider- 
able recent work on methods to engineer these states in 
optical modes [13, H, HI] . 

One can prepare a noon state with arbitrary N in the 
following way. We first prepare each oscillator in a super- 
position of two coherent states. Such a superposition is 
usually referred to as a Schrodinger cat state. One of the 
coherent states is the vacuum, and the other, which we 
denote |a), we choose so that its number-state distribu- 
tion is peaked near N phonons. Cat states are generated 
by applying the non- linear Hamiltonian (a^a) 2 , and this 
could be done for a nano-resonator using the techniques 
in [33| . Note that so long as JV > 1 , the probability that 
each oscillator is the zero-phonon state is much greater 
than the probability that it is in the one-photon state. 
We now make a measurement of the combined phonon 
number of the two oscillators. If the result of the joint 
energy measurement is N phonons, then, up to a nor- 
malization constant, the resonators are projected into the 
combined state 

\tp) oc c c N \N) + cic N -i\l) N + c 2 c N - 2 \'2,)n 

where we have used the fact that d n = c n . Since 
CN-m < cn for 1 < m < N, and in particular c rn <C cq 
for m < N/2, the state is a noon state with high fi- 
delity. This fidelity increases with N. If we wish to fix 
the value of N, we can decrease c%, and thus increase 
the fidelity, by squeezing the non-vacuum coherent com- 
ponent. If we choose a = s/~N, then the probability 
with which we obtain the desired measurement result is 
P(N) « (2-kN)- 1 / 2 . While this does decrease with N, 
it is quite reasonable for large values of N: for N = 50 
the measurement succeeds 1 time in 17 tries. In addi- 
tion, if one merely wants to produce a noon state with 
N > -/V min , then one can achieve a near- unity success rate 
by choosing a » V^min- For the purposes of graphically 
representing the noon states we define a subsystem of 
the two oscillators whose basis states are indexed by the 
phonon-number difference between the two oscillators. 
This basis is the \m) n, with m = 0, . . . , N. We can then 
plot the Wigner function for the state of this subsystem 
when the oscillators are in the noon state. In Fig. [3]^a) 
we plot this Wigner function for the noon state generated 
using the above procedure when a = 6 and the measure- 
ment projects onto the subspace with 20 phonons. 

We now turn to the task of preparing phase states. 
The problem of defining a phase operator conjugate 
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FIG. 3: Wigner functions for (a) a "noon" state with N = 20 
and (c) a phase-state with N = 20, generated using a QND 
measurement, (c) Inset: phonon number-difference distribu- 
tion for the same state. For these plots luminosity is pro- 
portional to the real part of the Wigner function, (b) Typi- 
cal phonon-difference distribution generated when the initial 
states are coherent. 



and 
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to the number operator goes back to Dirac [34[, 
was finally resolved in 1988 by Pegg and Barnett 
The phase eigenstates of the Pegg-Barnett phase oper- 
ator (the canonical phase states) are given by \6) oc 



limjv^oo J2o e 
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<G [0,27r]. For a finite max- 
imum phonon number N , there are N + 1 orthogonal 
phase states, \0 n ), corresponding to the phase angles 
8 n = 2irn/N . That is, when one is limited to a maximum 
excitation energy of fkuN, one can obtain a phase resolu- 
tion of A6 — 2n/N . The zero- phase state for iV-phonons 
is thus a real uniform superposition of all the number 
states from n = to N, and all the other 7V-phonons 
phase states are obtained from this by free evolution. 

Phase is not well-defined without a phase reference. 
While this statement is obvious for phase, it is, in fact, 
true for all physical observables, as was pointed out in the 
quantum context initially by Aharonov and Susskind [35| . 
For any observable to be physically meaningful, one must 
have a second system, referred to as a reference system, 
or reference frame, with which one can compare the sys- 
tem containing the observable. In the last few years 
the subject of reference frames for quantum systems has 
been studied extensively [3, [3(| . If one has two oscil- 
lators, then the state-space of these two systems "sup- 
ports" as a subsystem an oscillator with a well-defined 
phase. The number states of this oscillator are pre- 
cisely the number-difference states |m)jv defined above, 
and the zero-phase phase state of this oscillator is thus 
|#o) n = ELo \™)n/VN + 1. Along with reference sys- 
tems come the notion of a perfect reference state. This is 
a joint state of two systems, where the first is the base ref- 
erence system, such that the second system can be used 
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as a perfect (surrogate) reference for measuring a third 
system [lj|. The canonical phase state defined above, 
\8 ) N , is a perfect phase reference. It is also worth not- 
ing that of all the states spanned by the \rn)N, it is the 
canonical phase states that realize the highest entangle- 
ment between the two resonators. Thus the problem of 
generating phase states is also that of maximizing the 
entanglement generated by a joint energy measurement. 

The difficulty of preparing phase states is due to the 
fact that their perfectly flat number distribution is trun- 
cated sharply at a fixed number TV. However, the way 
in which a joint energy measurement folds together the 
states of the individual oscillators provides the ability 
to generate a flat distribution, and the subspace projec- 
tion allows one to generate the sharp cut-off. To begin 
we first recall how the joint energy measurement com- 
bines the coefficients c„ and d n of the states of the two 
resonators (Eq.([2])). This tells us that we need to choose 
c n dN-n = const., or if both oscillators have the same ini- 
tial state, c n CN- n = const, for some TV. This is satisfied 
by setting c„ oc e~^ n for some constant (3. We cannot ar- 
range for this to be true exactly, but we can make it true 
to a good approximation by choosing the initial states of 
the oscillators to be a superposition of a small number of 
squeezed coherent states. As examples, it turns out that 
high accuracies can be achieved for TV = 10 and TV = 20 
phase states with superpositions of only two and three 
squeezed states, respectively. Superpositions of multiple 
squeezed states can once again be obtained by applying 
the Hamiltonian (a^a) 2 for specified times. (We choose 
all the states to have the same squeezing, as this simplifies 
preparation (33|). The more coherent states in the initial 
superposition, then the higher the fidelity with which we 
can generate a phase states for a given TV. 

Defining the error in the preparation as / = 1 — 
|(V'I^o)n| 2 , where is the state we prepare, we find 
that using a superposition of two mildly squeezed states 
we can prepare an TV = 10 phase state with an error of 
less that 10~ 5 , and with three squeezed states we can 
prepare an TV = 20 phase state with / = 1.2 x 10~ 5 . To 
obtain these results we performed a numerical optimiza- 
tion over initial states. For generating an TV = 10 phase 
state the optimal amplitudes of the squeezed components 
are a.\ = 1.162, 012 — 3.277 with a (very mild) squeezing 
parameter of s = —0.097. For an TV = 20 phase state 
these are a± = 1.241, eti = 3.100, a$ — 5.024 with a mild 
squeezing parameter s = —0.1131. For higher values of 
TV one requires more coherent components in the initial 
state in order to maintain the same fidelity. 

As a final note, it is clear from Eq.([2]) that joint-energy 
QND measurements are universal for preparing states of 
the virtual oscillator (entangled states whose Schmidt ba- 
sis is {|^)at}) when combined with a universal technique 
for preparing states of a single oscillator. That of Law 
and Ebcrly [201 ] would be suitable for this purpose when 
TV is not too large. 
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